Vortex-Model for the Inverse Cascade of 2D-Turbulence 
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We generalize Kirchhoff's point vortex model of two-dimensional fluid motion to a rotor model 
which exhibits an inverse cascade by the formation of rotor clusters. A rotor is composed of two 
vortices with like-signed circulations glued together by a nonelastic bond. The model is motivated 
by a treatment of the vorticity equation representing the vorticity field as a superposition of vortices 
with elliptic Gaussian shapes of variable widths, augmented by a suitable forcing mechanism. The 
rotor model opens up the way to discuss the energy transport in the inverse cascade on the basis of 
dynamical systems theory. 
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The theoretical treatment of the longstanding problem 
of turbulent flows [l|-[H has to relate dynamical systems 
theory with nonequilibrium statistical physics [5]. The 
central notion of physical turbulence theory is the con- 
cept of the energy cascade, highlighting the fact that tur- 
bulent flows are essentially transport processes of quan- 
tities like energy or enstrophy in scale. Although well- 
established theories due to Richardson, Kolmogorov, On- 
sager, Heisenberg and others (for reviews we refer the 
reader to (HfS]) can capture gross features of the cas- 
cade process in a phenomenological way, the dynamical 
aspects are by far less understood, and usually are in- 
vestigated by direct numerical simulations of the basic 
fluid dynamic equations. An exception, in some sense, 
are inviscid fluid flows in two dimensions. Based on the 
work of Hclmholtz [6] Kirchhoff [7] pointed out that the 
partial differential equation can be reduced to a Hamilto- 
nian system for the locations of point vortices, provided 
one considers initial conditions where the vorticity is a 
superposition of delta-distributions (we refer the reader 
to the works of Aref [las well as the monographs 
[icil|). Due to Onsager [lTjj (for a discussion we refer the 
reader to [12), a statistical treatment of point vortex dy- 
namics is possible for equilibrium situations because of 
the Hamiltonian character of the dynamics, provided the 
ergodic hypothesis holds. Extensions to nonequilibrium 
situations based on kinetic equations have been pursued 
by Joyce and Montgomery [13| . Lundgren and Pointin 
14 1 , as well as more recently by Chavanis [l5[ . 

The purpose of the present Letter is to generalize 
Kirchhoff's point vortex dynamics to a rotor model which 
exhibits an inverse cascade, i.e. the spontaneous forma- 
tion of large scale velocity fields from an initially random 
distribution of vortices. Kraichnan [16| has emphasized 
that such an inverse cascade can exist in two-dimensional 
flows. The predicted scaling of the energy spectrum has 
been experimentally verified by Paret and Tabel ing [l7| . 
Recently, it has been pointed out by Chen et al. [18[ that 
a vortex thinning mechanism is an important feature of 
the inverse cascade. A similar mechanism is contained in 
the rotor model presented in this Letter. We will explic- 
itly demonstrate that the large scale fields are generated 



by a clustering of rotors, induced by a change of the ro- 
tor size due to a large scale straining field, i.e. a kind 
of vortex thinning process. In a second part of the Let- 
ter we shall discuss how the rotor model is related to 
the vorticity equation of dissipative fluid dynamics with 
forcing. 

The Model. Instead of single vortices located at posi- 
tion Hi with circulations Ti we consider vortex couples, 
where each vortex of the couple i has circulation F^, lo- 
cated at the positions Xj and yj. In the following we 
denote these objects as rotors. We assume that the posi- 
tions xj and yi are determined by the evolution equations 



Xj = -(D - |xj - y l |)e i + r j u(x i - yj) 

+ r j"t u ( Xi _ x i) + u(x * _ Yj)] 

j 
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ji = (A - |y« - Xi|)ej + T i u(y l - Xi) 

+ Yl T ji u (Yi - Yi) + "(y* - x,)] (1) 



We have defined the unit vector e ? - 



The ve- 



|x,-y,r 

locity field u(r) is the velocity field of a point vortex 
centered at the origin, u(r) = e z x For vortices 

moving inside a closed regime, the velocity field has to 
be changed based on the introduction of mirror vortices, 
as is described e.g. in [§], The size of a single rotor 
relaxes with relaxation time 2/j to Dq. 

In our rotor model, the vortex couple is glued together 
by an inelastic coupling. The rotors generate a far field 
similar to the one of an elliptical vortex. Furthermore, 
due to the inelastic coupling it possesses an internal de- 
gree of freedom, which is sensitive with respect to a shear 
velocity field, as is easily seen from the multipole expan- 
sion of the evolution equation ((T|) : 



Ti = 7(A) - n)- + 2r i u(r i ) + T 3 r * ■ Vu(Ry) (2) 
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R* = ^21X11^) (3) 
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Here, we have introduced the center R, 



of a 



rotor, the distance vector R^- = Ri — R^ between two 
rotors, as well as the distance vector Yi = Xj —yt, which is 
a measure of the elliptical character of the velocity field of 
the rotor. We have included only the leading terms in the 
multipole expansion of the velocity field with respect to 
|r|/|R|. As is seen from eq. ([2]) the elliptical deformation 
of the velocity field, i.e. the size and the direction of 
Ti, depends on the shear induced by the other rotors. 
The inelastic coupling between the vortices of the couple 
can be considered to model viscous damping as well as 
forcing. 
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FIG. 1: Inverse cascade in the rotor model with TV = 200 
rotors : a) An inital stage starting from random initial con- 
ditions. 

b) Evolved state showing the formation of rotor clusters. 

Numerical results We have numerically solved the dy- 
namical system (|T|). The temporal evolution of 200 rotors 
with an equal number of positive and negative circula- 
tions starting from a random initial condition exhibits 
an inverse cascade. Large scale vortical motion is gen- 
erated by a clustering of rotors. A typical time series is 
exhibited in figs. 1, a,b , for the parameter values (7 = 5., 
Dq = .5, T = ±1, R = 20). The temporal evolution of a 
number of rotors with the same circulations is exhibited 
in figs. 2, a,b. Starting from random initial positions of 
the rotors a fluctuating lattice of rotor clusters appear. 

The fact that the rotor vortex system exhibits a pro- 
nounced inverse cascade already for moderate numbers 
of rotors (the figures contain 100 rotors) allows us to in- 
vestigate the inverse cascade using methods of nonlinear 
dynamics. 

Furthermore, the rotor model incorporates the aspect 
of vortex thinning, due to a possible change of the ellip- 
ticity of the rotor in much the same way as identified in 
the experiments of Chen et al. [l8| . Hence, it is a mini- 
mal dynamical model containing the mechanisms of the 
inverse cascade. In the following we shall discuss the ori- 
gin of the formation of clusters of rotors with like-signed 
circulations. 
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FIG. 2: Inverse cascade in the rotor model (N=200) with 
like signed circulations: a) Inital stage starting from random 
initial conditions, b) Evolved state showing the formation of 
rotor clusters generating large scale flow. 



Interaction of two rotors. It is straightforward to show 
that the center of vorticity Ti^i+r^) is a conserved 
quantity. The distance vector R between the two rotors 
obeys the evolution equation ([3]), written explicitly 
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For the following it is convenient to introduce the ab- 
breviation 



h(ei,Tn) = e* • [e z x e fl J(e !: • e R ) 



(5) 



with the unit vectors = e# = We obtain the 
equation for the relative distance 



R 



{ Il±lA^[rlh{e u e r )+rlh{ e2ier )] (6) 



The evolution equation for the relative distance of a 
rotor reads 



ri 



7(A) - n) h 2— e, x 
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We have to determine the quantities r 2 , r 2 , which are 
determined by the evolution equations 



ri =j(D -ri)-2-^-^h(e 1 ,e R ) 



(8) 



We can solve iteratively for small deviations of rj from 
D : 



n = D -2 



£2 Do 
2tt R 2 



e -T( t - i ')/ l (e 1 (i'),e fl (t')) (9) 



A similar treatment applies to r 2 . Assuming that the 
damping constant 7 is large compared to the rotation 
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frequency of the rotor, adiabatic approximation of the 
integral yields to lowest order in 7 _1 

r? = £>g[l- 4^^(61, eji )] (10) 

and a similar equation for r\. Here, the last terms on 
the right hand side arise due to the change of the size 
or the rotors, connected with a change of the far field, 
induced by the mutually generated shear. It mimics the 
mechanism of vortex thinning, identified in [Lsj . 

The relative motion of the rotors obeys the evolution 
equation 



R = 



(ri + r 2 ) d, 

2ir R 3 

a r 2 



°{[h(ei,r R ) + h(e 2 ,T R ) 
„ r i 



2-k^R 2 



h(ei,r R 



2irjR 2 



h(e 2 ,r R y] (11) 



We now average the evolution equation with respect 
to the rotations of the vectors ej(i), i = 1,2, taking into 
account that the average (h(ei,r R )) = vanishes. Fur- 
thermore, the average (h(ei,r R ) 2 ) = a is positive. As a 
consequence, the relative distance behaves according to 



R 



(r 1 + r 2 ) 2 Dl 

" (2tt) 2 7 i? 5 ' 



(12) 



Two rotors approach each other, except for T\ = — r 2 . It 
is important to stress that this attractive relative motion 
arises only if we include the irreversible effect of the strain 
induced stretching of the rotors. 

Relation to vorticity equation. In the following we shall 
discuss the relationship of our point vortex model to the 
vorticity equation of fluid dynamics. Thereby, we replace 
the point vortices by vortices with elliptic Gaussian pro- 
files, along similar lines as described by Mclandcr, Sty- 
czek, and Zabusky [l9j . 

We consider the two dimensional vorticity equation in 
Fourier-space 

u(k,t)-ik-J rfk'u(k>(k-k',t)w(k' ) t) = -vk 2 u(k,t) 

(13) 

with u(k) = -j^y [e z x py]. In order to deal with vortex 
dynamics, we perform the ansatz 



; ik-X;j(t)-(-W}(k,t) 



(14) 



For Wj(k,t) = we obtain the field of point vortices 
w(x, t) — Tj5(x — Xj(t)). We now insert our ansatz 
(1141) into the vorticity equation and obtain 



3 ik-Xj(t)+W 3 -(k,t) 



ik-xj(t) +Wj(k,t) + vk 2 



ik ■ r 3 T i J dk'u(k')e 4(k - k ' ) x J +tk ' Xi 

,w>(k-k',*)+wi(k',t) / 15 ^ 
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This equation contains sweeping of the vortices, encoded 
in the location of the centers Xj (t) , as well as a change 
of the shapes of the vortices, contained in Wj(k,t). 

A separation of these two processes becomes possible 
by defining sweeping via the terms in the evolution equa- 
tion (|15[) which are proportional to k: 



(16) 



We have defined the velocity fields 

Ujj(r) = J dk'u(k')e- lk ' r e w ^- k '^ +Wl< - k '^ (17) 

This is the extension of the set of evolution equations for 
the (5-point vortices. Let us now consider the resulting 
equation for the shapes, which read 



Wi(k,t) = -vk 2 

+ik-J2 T i J rfk'u(k')e- lk ' [: 



(18) 



-x,] e W i (-k',t)+Wi(k',i) 



x e 



Wi(U-W ,i)-Wi(-k' ,t)-Wi(k,t) _ i 



Here, the sum includes also the terms with i 
now invoke the approximation 



Wi(k,t) 



~kd(t)k 



I. We 



(19) 



introducing the symmetric matrix d. Approximating 
the last term of the righthand side of eq. (fT5|) by 
i[kCik' + k'Qk] we arrive at 

Ci = vE + TiiSuCi + CiST] (20) 
+ ^T^Suixi - x t )Ci + CiSuixi - X/ ) T ] 

Here, we explicitly have introduced the matrix Su = 
[Wu(xi — x;)] and have singled out the term with I = i. 

In order to make contact with the previously intro- 
duced rotor model we determine the velocity field eq. 
(fTT]) up to the first order in d + Ci valid for widely sep- 
arated vortices 

dk' ■ u(k')e' lk '- r e-2k'(c J +c 1 )k' 

« e 2 x[l + iv r (Q+Q)V r ]^ 

The evolution equation for the vortex centers then read 

Xi - x ; 



Xj - X; 

47r|x l - xA 



(21) 



x e 



This equation is identical to the equation ([3]) for the vor- 
tex centers of our point vortex model, provided we rep- 
resent the matrix as the direct product Ci = rj(£)rj(t). 
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This corresponds to an infinitely thin elliptical vortex, 
oriented in redirection. In this case, the evolution equa- 
tion for Ti, obtained by inserting Gj = i"j(t)rj(i) into eq. 
([2"U|). is identical to the model equation eq. © of the 
rotor model with 7 = 0. 

As we have seen the forcing term 7 (.Do — fi)^i is cru- 
cial for the emergence of the inverse cascade in the rotor 
model (|T|). It is tempting to include forcing and damping 
in the following way: 

T t = -aTi + fi 

x, = ]rr ; U. i( (x 4 -x ( ) + U^) (22) 
I 

di = vE + ^Co-Cij + riiSuCi + CiST:] 

+ y^^Ti[Su(xi -xi)Ci +CiSu(xi -x ; ) T ] 
i 

Such type of forcing may be obtained from the vortic- 
ity equation (| 13[) by just adding a linear damping term, 
— aoj(k, t) as well as the forcing term F(k, t), 



(23) 



x [f j + iTjk ■ Uj (t) - -r,k 7 (Co (t) - C, (t)) k] 

where the centers Xj = Xj(i) + ^Uj(i) as well as the 

shapes Wj = Wj — ^^-jk(C (t) — Cj(t))k are close to 
the centers and the shapes of the elliptical vortices. The 
first contribution in eq. (|2"3")l leads to a modulation of the 
circulation, the second term describes a shift of the rotor 
center and the third one corresponds to a modification 
of the width of the Gaussian vortex shape counteracting 
the viscous broadening of the vortex. 



Conclusions We have presented a generalized point 
vortex model, a rotor model, exhibiting an inverse cas- 
cade based on clustering of rotors. We have discussed 
how this rotor model can be derived from the vorticity 
equation by an expansion of the vorticity field into a set of 
elliptical vortices at locations Xj(£) and shapes C,(t). An 
important point has been the inclusion of a forcing term, 
which prevents the elliptical far field of the rotors from 
diffusing away. The added forcing term breaks the sym- 
metry Ti — > —Ti, t — > —t. This symmetry breaking lies at 
the origin of cluster formation and the inverse cascade, as 
can be seen from the two-rotor interaction inducing in av- 

D 2 

erage a relative motion proportional to — + T 2 ) 2 / R 5 . 

The presented rotor model can be investigated by 
applying methods from dynamical systems theory like 
the evaluation of finite time Ljapunov exponents and 
Ljapunov vectors. The model system (|2"2"|) may also be 
studied as a stochastic system by considering the veloc- 
ity Ui(i) to be a white noise force. The correspond- 
ing Fokker-Planck equation allows one to draw analogies 
with quantum mechanical many body problems. Fur- 
thermore, we emphasize that a continuum version of the 
model equations (|22j) leads to a subgrid model exhibiting 
analogies with the work of Eyink [20j . 

It will be a task for the future to investigate the clus- 
ter formation from a statistical point of view, based on 
the formulation of kinetic equations, along the lines as 
has been performed for fully dev elop ed turbulence [2lj |. 
and Raylcigh-Benard convection [22J . In this respect we 
hope to find a relation to the kinetic equation for the 
two-point vorticity statistics recently derived on the basis 
of the Monin-Lundgren-Novikov hierarchy, taking condi- 
tional averages from direct numerical simulations 23 1. 
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